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Cirsium arvense populations showed 2n = 2x = 34 chromosome number, forming 
mostly bivalents, although some amount of univalents and quadrivalents were also 
formed. The Ghasemloo populations showed the highest value of total and terminal 
chiasmata, while the Khalkhal population showed the least values of the same parameters. 
The Meyandoab population had the highest value of intercalary chiasmata. Meiotic 
abnormalities including, laggard chromosomes, stickiness, micronuclei formation, 
unorganized chromosomes, multipolar cell formation and desynapsis were observed 
almost in all the populations studied. The occurrence of large pollen grains (possibly 2n 
pollen grains) was observed along with smaller (normal) pollen grains in some of the 
populations. B-chromosome (0-1) was observed only in the Taleghan population. RAPD 
and ISSR primers used produced many reproducible bands. Some of the populations 
showed the presence of specific bands, while some bands were present in all except one 
population. UPGMA and NJ trees of RAPD and ISSR data grouped the populations in two 
different clusters/groups, indicating their genetic difference which is discussed in details. 
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The genus Cirsium Mill. {Asteraceae) 
contains about 250 perennial, biennial or rarely 
annual and spiny species (Zomlefer 1994, 
Bures et al. 2004). Growing mainly in Northern 
Hemisphere (Europe, North Africa, Siberia, 
Central and East Asia, West, East and Central 
America). 

The genus is considered to be complex in 
nature partly due to high degree of hybridization 
occurring among different species and also 
because of highly variable germplasm interacting 
with different environmental conditions. In 
fact several hybrids have been described from 


Caucasus and the adjoining regions of Asia 
Minor (Bures et al. 2004). About 44 Cirsium 
species have been reported in Flora Iranica 
(Rechinger 1979), which have been classified 
in five sections. Hybrid plants persist primarily 
through vegetative growth, forming clusters of 
flowering shoots often connected by rhizomes. 
Because hybrids are usually fertile, they can 
often produce introgressive hybrids with the 
parental species or triple hybrids with other taxa 
(Bures et al. 2004). 

Limited cytological studies have been 
performed on the genus Cirsium (Frankton 
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Fig. 1. Map of Iran, showing localities of the Cirsium arvense populations studied. Population code. 1. 
Asara Karaj. 2. Porkan Karaj. 3. Kordan Tehran. 4. Taleghan Ghazvin. 5. Malard Karaj. 6. Shahrestanak 
Karaj. 7. Tafresh Markazi. 8. Mahallat Markazi. 9. Kashan. 10. Ghasemloo West Azarbayejan. 11. 
Tabriz West Azarbayejan. 12. Klialklial Ardebil. 13. Meyandoab West Azarbayejan. 



and Moore 1961, Tonian 1981a, 1981b, Tonian 
1982, Dempsey 1994, Ghaffari 1999, Ozcan 
et al. 2008) and only recently some cytological 
study concerned with chromosome pairing 
and B-chromosomes at the species level 
have been reported from Iran (Nouroozi et al. 
2010, 2011). In these studies, some indication 
of the occurrence of unreduced (2n) pollen 
grains in other Cirsium species but not in C. 
an’ense has been provided. Since C. arvense 
is of wide distribution in Iran and grows in 
various environmental conditions, we planned 
to study some of these populations for the 
potential of unreduced pollen grain formation 
and B-chromosome occurrence. Moreover the 
present study reports molecular diversity of 13 
populations of C. an’ense growing in Iran for the 
first time. 


Different molecular markers including 
RAPD (Random Amplified Polymorphic 
DNA) as well as ISSR (Intersequence Simple 
Repeats) have been used for studying genetic 
diversity and identification of hybrids different 
plant species (e.g., Chen et al. 2004, Sheidai et 
al. 2007, Wei et al. 2008, Tafvizei et al. 2010, 
Sheidaei et al. 2010), including Asteraceae taxa 
(Ash et al. 2003, Dogan et al. 2007). 

Materials and Methods 
Plant materials 

Thirteen populations of Circium arvense 
(L.) Scop, were collected for cytological and 
molecular studies (Fig. 1). 

Cytological studies 

Meiotic studies were performed on young 
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flower buds collected using minimum 100 
metaphase/diakinesis pollen mother cells 
(PMCs) and 500 anaphase and telophase cells 
for data collection (Sheidai and Rashid 2007). 
Pollen stainability as a measure of fertility was 
determined by staining minimum 1000 pollen 
grains with 2% acetocarmine: 50% glycerin 
(1:1) for about 1/2 hr. Round complete pollens 
which were stained were taken as fertile, while 
incomplete, shrunken pollens with no stain 
were considered as infertile (Sheidai and Rashid 
2007). 

RAPD analysis 

Thirty decamer RAPD primers of Operon 
Technology (Alameda, Canada) belonging to 
OPA, OPH sets were used in this study. DNA 
extraction was done by using the CTAB method 
(Murry and Tompson 1980) with modification 
described by De La Rosa et al. (2002). The PCR 
reaction mixture consisted of 1 ng template 
DNA, 1 x PCR buffer (10 mM Tris-HCL pH 
8.8, 250 mM KCL), 200 pM dNTPs, 0.80 
pM 10-base random primers and 1 unit of Taq 
polymerase, in a total volume of 25 pL. DNA 
amplification was performed on a palm cycler 
GP-001 (Corbet, Australia). Template DNA was 
initially denatured at 92°C for 3 min, followed 
by 35 cycles of PCR amplification under the 
following parameters: denaturation for 1 min 
at 920C, primer annealing for 1 min at 36°C 


and primer extension for 2 min at 72°C. A final 
incubation for 10 min at 72°C was performed 
to ensure that the primer extension reaction 
proceeded to completion. 

The PCR amplified products were separated 
by electrophoresis on a 2% agarose gels using 
0.5 x TBE buffer (44.5 Mm Tris/Borate, 0.5 
Mm EDTA, pH 8.0) or 6% polyacrylamide gels. 
The gels were stained with ethidium bromide 
and visualized under UV light (Sambrook et 
al. 2001). A 100 bp DNA ladder (GeneRuler, 
Fermentas) was used as the molecular standard 
in order to confirm the appropriate RAPD 
markers. These markers were named by primer 
origin, followed with the primer number and the 
size of amplified products in base pairs. 

ISSR assay 

Total genomic DNA was extracted from 
fresh leaves using the CTAB method by Murry 
and Tompson (1980) with the modification 
described by De la Rosa et al. (2002). Six 
ISSR primers used are (GA) 9 T, UBC810, 
UBC811, UBC834, UBC849 and CA7GT 
commercialized by UBC (the University 
of British Columbia). PCR reactions were 
performed in a 25 pL volume containing 10 
mM Tris-HCl buffer at pH 8; 50 mM KCl; 1.5 
mM MgCl 2 ; 0.2 mM of each dNTP; 0.2 pM of 
a single primer; 20 ng genomic DNA and 3 unit 
of Taq DNA polymerase (Bioron, Germany). 


Table 1. Cytogenetic characters in the Cirsium arvense populations studied 


Locality 

2n 

ROD 

RB 

I 

IV 

IX 

TX 

TOX 

B 

Taleghan 

34 

2.08 

7.80 

5.50 

0.80 

0.88 

20.86 

21.75 

0-1 

Tabriz 

34 

3.30 

5.90 

4.90 

1.38 

1.9 

20.20 

22.10 

0 

Ghasemloo 

34 

3.80 

9.50 

2.60 

0.50 

1.8 

24.87 

26.60 

0 

Tafresh 

34 

3.93 

8.76 

2.26 

0.86 

2.23 

24.40 

26.56 

0 

Mahallat 

34 

2.32 

11.24 

3.44 

0.04 

1.44 

24.96 

26.40 

0 

Khalkhal 

34 

4.23 

6.73 

5.80 

0.13 

1.60 

18.23 

19.90 

0 

Meyandoab 

34 

3.70 

9.10 

3.60 

0.28 

2.60 

22.60 

25.35 

0 


Abbreviations. ROD. Mean of rod bivalents. RB. Mean of ring bivalents. I. Mean of univalents. IV. Mean of 
quadrivalents. IX. Mean of intercalary chiasmata. TX. Mean of terminal chiasmata. TOX. Mean of total chiasmata. B. 
B-chromosomes. 




Fig. 2. Representative meiotic cells in the Cirsium arvense populations. A-C. Metaphase cells in the Tabriz, Tafresh 
and Meyandoab populations. D, E. Metaphase cells showing chromosomes stickiness in the Khalkhal and Tabriz 
populations. F, G. Meiocytes showing complete desynapsis in the Khalkhal and Tafresh populations. H-J. Multipolar 
cells in the Tafresh, Khalkhal and Mahallat populations. K, L. Potential unreduced pollen grains (larger size pollen 
grains) in the Khalkhal and Meyandoab populations. Scale bar= 10 pm. 


Amplifications reactions were performed in 
Techne thermocycler (Germany) with following 
program: 5 min initial denaturation step 94°C, 
30 s at 94°C; 1 min at 50°C, 1 min at 72°C. 
The reaction was completed by final extension 
step of 7 min at 72°C. Amplification products 
were visualized by running on 2% agarose gel, 
following ethidium bromide staining. Fragment 
size was estimated by using a 100 base pairs (bp) 
molecular size ladder (Fermentas, Germany). 


Data analysis 

In order to determine significant difference 
in meiotic characteristics, paired sample y} 
test was performed. Principal components 
analysis (PCA) was performed to show meiotic 
differences among the populations studied. 

RAPD and ISSR bands obtained were treated 
as binary characters and coded accordingly 
(presence =1, absence = 0). Jaccard similarity 
as well as Nei’s genetic distance (Nei 1973) 
was determined among populations studied 
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and used for grouping of the genotypes by 
UPGMA (Unweighted Paired Group with 
Arithmetic Average) and NJ (Neighbor-Joining) 
clustering methods and ordination based on 
principal coordinate analysis (PCO) and PCA 
(Podani 2000, Weising et al. 2005). NTSYS 
Ver. 2.02 (1998) and DARwin ver. 5 (2008) 
was used for clustering and PCO analyses. In 
order to determine molecular difference between 
populations, AMOVA test was performed. For 
test of agreement between RAPD and 1SSR trees 
obtained. Mantel test was performed (Podani 
2000) by GENALEX 6 (Peakall and Smouse 
2006). 

Results 

Cytology 

Cirsium arven.se populations showed 2n = 
2x = 34 chromosome number (Fig. 2, A-C), 
supporting Nouroozi et al. (2010, 2011). All 
populations formed mostly bivalents, although 
some amount of univalents and quadrivalents 
were also formed (Table 1). 

The Ghasemloo populations showed the 
highest value of total and terminal chiasmata 
(26.60 and 24.87, respectively), while the 
Khalkhal population showed the least values 
of the same parameters (19.90 and 18.23, 
respectively). The Meyandoab population had 
the highest value of intercalary chiasmata (2.60). 

Various meiotic abnormalities including, 
laggard chromosomes, stickiness, micronuclei 
formation, unorganized chromosomes, 
multipolar cell formation and desynapsis were 
observed almost in all populations studied (Table 
1, Fig. 2, D-L). 

Chromosomes stickiness occurred in both 
metaphase of meiosis 1 and II. This occurred 
between two bivalents or more, forming 
complete clump in some of the cases (Fig. 1, 
D, E). The highest percentage of metaphase 
I cells showing chromosome stickiness was 
observed in Tabriz population (42%), followed 
by Meyandoab population (25%). Taleghan 
population was the only one not showing 


metaphase I stickiness. 

Similarly the Mahallat populations showed 
the highest percentage of metaphase II cells 
showing chromosome stickiness (20%), 
followed by the Taleghan population (15%). The 
Tabriz population was the only one not showing 
metaphase II stickiness. Complete chromosome 
clumping was mostly observed in the Taleghan, 
Ghasemloo and Taffesh populations. 

Different degrees of univalent formation 
to complete desynapsis (all univalents) were 
observed in the populations studied (Fig. 2, F, 
G). Chromosomes stickiness during anaphase 
I and II leading to the formation of anaphase 
bridges in some of the cases were observed 
in. The number and thickness of these bridges 
varied among different meiocytes. Laggard 
chromosomes of different numbers occurred 
in Taleghan and Mahallat populations. These 
laggards may be the reason for micronuclei 
formation in telophase I and II. 

Anaphase abnormalities leading to the 
formation of multipolar cells were observed 
in all populations (Fig. 1, H-J). The highest 
percentage of triads occurred in the Ghase ml oo 
population followed by the Taleghan population. 
This latter population also showed the highest 
percentage of pentads and hexads. % 2 test showed 
significant difference for meiotic abnormalities 
among the populations studied. 

The occurrence of large pollen grains 
(possibly 2n pollen grains) was observed along 
with smaller (normal) pollen grains in five 
populations of Taleghan, Tabriz, Ghasemloo, 
Tafresh and Meyandoab (Fig. 1, K, L), all 
showing triad formation (Fig. 1, H, J). The large 
pollen grains comprised about 1% of pollen 
grains in these populations. B-chromosome (0- 
1) was observed only in the Taleghan population. 

Chi square test showed that the Kha lk hal 
population has significantly lower chiasma 
frequency compared to that of the Tafresh, 
Mahallat and Ghasemloo populations (yj = 
14.70, p<0.01). Significant difference was 
observed for meiotic abnormalities among 
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Fig. 3. PCAplot ofmeiotic data in the Cirsium arvense populations. 


populations studied too. 

PCA plot showed cytogenetic differences 
of the populations studied (Fig. 3). Two 
populations of Taleghan and Mahallat are placed 
in different positions in the PCA plot, far from 
the other populations. Similarly Tabriz and 
Khalkhal are placed far from the others. 

RAPD analysis 

Out of 20 RAPD primers used, 13 primers 
produced reproducible bands (Fig. 4). In total 
178 RAPD bands (loci) were obtained out of 
which only 8 bands ranging in size from 500- 
2000 bp were common in all 13 populations 
studied and 180 bands were polymorph. 

Among the primers used OPC-08 produced 
the highest number of bands (23), while primers 
OPI-12 produced the lowest number of bands 
(10). The Primer OPC-08 also produced the 
highest number of polymorphic bands (23). 


The primer OPC-08 produced 6 specific 
bands, primer OPC-10 produced 5 specific 
bands, while some other primers produced 2 and 
3 specific bands and the primers OPR-08, OPI-1 
and OPC-05 produced 1 specific band. 

Some of the populations showed the presence 
of specific bands, for example only Kordan 
population had band 1400 bp of RAPD primer 
OPI-12, band 1100 bp of the primer OPC-05 and 
band 1300 bp of the primer OPC-10. 

The Tabriz population showed specific bands 
of 400 bp of the primer OPR-08, while Asara 
population was the only population having band 
700 bp of the primer OPM-10, band 900 bp 
of the primer OPC-08 and 300 bp of OPA-02. 
Kashan population had the band of 500 bp of the 
primer OPA-02. 

Some bands were present in all except one 
population, for example bands 690 bp of the 
primer OPM-10 and 550 bp of OPC-05 were 
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Fig. 4. OPA04 RAPD primer profile of the Cirsium arvense populations. Abbreviations. L. Ladder. 1. Asara. 2. 
Porkan. 3. Kordan. 4. Taleghan. 5. Malard. 6. Shahrestanak. 7. Tafresh. 8. Ghasemloo. 9. Mahallat. 10. Tabriz. 11. 
Kashan. 12. Klialklial. 13. Meyandoab. N. No DNA. 


Tabriz Ghasemloo 


Meyandoab 


Shahrestanak 


Taleghan Kordan 


Fig. 5. NJ tree of RAPD data in the Cirsium arvense populations. 


absent only in Taleghan population, the bands 
550 bp of the primer OPC-05, bands 650 bp 
of the primers OPC-08 and OP 1-12 were only 
absent in Kashan population. 

UPGMA and NJ dendrograms of RAPD 
data produced similar results supported by PCO 
ordination plot (Figs. 5, 6). The Cophenetic 
correlation of NJ tree was higher (r = 0.98) and 


therefore it is discussed bellow. 

In general two major clusters were obtained. 
The populations of Asara, Shahrestanak, Porkan, 
Malard, Kordan and Taleghan formed the 
first major cluster, out of which, Kordan and 
Taleghan show more RAPD similarity and the 
Asara population differ the most from the other 
populations. 
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Fig. 6. PCAplot of RAPD data in the Cirsium arvense populations. 


Seven other populations form the second 
major cluster, including Taffesh, Ghasemloo, 
Tabriz, Meyandoab, Mahallat, Kashan and 
Kha lk h a l. Out of which, Kashan and Kha lk hal 
and also Ghasemloo and Tabriz show more 
similarity in RAPD while the Tafresh population 
differ the most from the other populations. 

PCA plot reveals RAPD molecular diversity 
of the populations in more details (Fig. 6). The 
Porkan population stands far from the other 
populations, to which the Asara population joins 
with some distance. Meyandoab stands far in the 
other side of the plot, to which the Khalk hal and 
Tabriz populations are joined. 


ISSR analysis 

The six ISSR primers used in this experiment 
produced 83 reproducible bands (Fig. 7), out 
of which 80 bands were polymorph and 3 
bands were monomoiph. The highest number 
of polymorphic bands was obtained for the 
ISSR primer CA7G (20), while ISSR primers 
UBC811 produced only 9 polymorphic bands. 

In total 16 specific bands were obtained, 
with ISSR primer (GA) 9 T producing the highest 
No. (5) and primer UBC811 producing only 1 
specific band. Some of the populations showed 
the presence of specific ISSR band/locus, for 
example, Shahrestanak population was the only 
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Fig. 7. CA7GT ISSR profile of the Cirsium arvense populations. Abbreviations: N. No DNA. 1. Meyandoab. 2. 
Khalkhal. 3. Kashan. 4. Tabriz. 5. Mahallat. 6. Ghasemloo. 7. Tafresh. 8. Shahrestanak. 9. Malard. 10. Taleghan. 
11. Kordan. 12. Porkan. 13. Asara. L. Ladder. 



Meyandoab 


Fig. 8. NJ tree of ISSR data in the Cirsium arvense populations. 


population having ISSR band 2100 bp of the primer UBC810. The Porkan population was 
primer UBC849, band 1500 bp of the primer the only one showing ISSR band 2300 bp of 
(GA) 9 T and bands 2000 and 5000 bp of the the primer UBC810, bands 760, 1100 and 200 
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Fig. 9. PCAplot of ISSRdata in the Cirsium aivense populations. 


bp of the primer CA7GT, band 110 bp of the 
primer (GA) 9 T. The Kordan population was the 
only population showing bands 2300 and 5000 
bp of the primer UBC810 and only Taleghan 
population had the band 600 bp of the same 
primer. 

Some 1SSR bands were present in all 
populations except one, for example, 1SSR band 
300 bp of the primer UBC810 was present in all 
populations except in the Malard and Kashan 
population did not have the band 600 bp of the 
primer (GA) 9 T, which was present in the other 
populations. Meyandoab population did not have 
the bands 700 and 850 bp of the primer CA7GT, 
while the Malard population did not show the 
band 300 bp of the primer UBC810 occurring in 


the other populations. 

UPGMA and NJ dendrograms of 1SSR data 
produced similar results supported by PCO 
ordination plot. The Cophenetic correlation 
of NJ tree was higher (r = 0.96), therefore it is 
discussed bellow (Figs. 8,9). 

In general two major clusters were obtained, 
separating Porkan, Asara, Kordan, Shabrestanak, 
Taleghan and Malard in one cluster and the 
other populations in the second major cluster. 
Kordan and Shahrestanak in the first cluster and 
Ghase ml oo and Mahallat in the second cluster, 
show more similarity in 1SSR characteristics. 

Mantel test performed between RAPD and 
1SSR trees showed t = 0.03, p = 0.37, showing 
agreement between the two trees obtained. 
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AMOVA test performed between populations 
placed in the two major clusters in both RAPD 
and 1SSR markers showed significant difference 
(p<0.002) for 1SSR with 14% of variance 
between two groups and 86% within populations 
in each group, and (p<0.001) for RAPD with 
21% of variance between the two groups and 
79% within populations. 

Discussion 

Cytology 

All the populations studied showed 2n = 2x 
= 34 chromosome indicating that no polyploid 
population we could encounter. The populations 
studied differed in meiotic behavior and chiasma 
frequency. These results indicate a change in 
meiotic control over chromosome pairing and 
chiasma formation during species diversification. 
Variation in chiasma frequency and localization 
is genetically controlled (Quicke 1993) and has 
been reported in several plant species as well 
as in crop plant varieties (Rees and Dale 1974, 
Rees and Jones 1977). 

All the populations studied although 
are diploid, shows that the occurrence of 
quadrivalents as a result of heterozygote 
translocation. Such heterozygote translocations 
lead to the formation of new linkage groups. 
Therefore both mechanisms of a change in 
chiasma frequency and the occurrence of 
heterozygote translocations in the C. arvense 
populations with the same chromosome number 
may be a mean for generating new forms of 
recombination influencing the variability within 
natural populations in an adaptive way (Rees 
and Dale 1974). 

Molecular analyses 

The presence of RAPD and 1SSR 
polymorphic bands in the populations studied 
indicates the presence of genetic polymorphism 
in these genotypes. Moreover the occurrence 
of specific bands/loci only in some of the 
populations illustrates the occurrence of unique 
insertion/deletion in DNA material of these 


genotypes. 

As stated before, there were bands which 
occurred in some of the populations but were 
absent in the others. Since, even single base 
change at the primer annealing site is manifested 
as appearance or disappearance of RAPD and 
1SSR bands, these bands may indicate the 
occurrence of genetic changes in the genome of 
the populations through the loss or rearrangement 
of some of their nucleotides. Chromosomal 
crossing over during meiosis may result in loss 
of primer attachment pair sites in the offspring 
leading to novel RAPD pattern in the offspring 
(Smith et al. 1996). We also observed significant 
change in chiasma frequency and distribution as 
well as heterozygote translocations in C. arvense 
populations which may be responsible for 
RAPD/ISSR diversity. 

Both RAPD and 1SSR trees obtained are in 
general agreement separating the populations 
studied in two groups, indicating their genetic 
difference, which is also partly supported by 
cytogenetic analysis. It is interesting to say that 
populations of Asara, Porkan, Kordan, Taleghan, 
Malard and Shahrestanak are placed in three 
neighbor provinces of Tehran, Ghazvin and 
Alborze, which form the first group/cluster in 
RAPD and ISSR analyses, while, populations 
of Tafresh and Mahallat (Markazi Province), 
Ghasemloo and Meyandoab (West Azarbayejan), 
Tabriz (East Azarbayejan, Khalkhal (Ardebil) 
and Kashan (Isfahan) which comprise the 
second group/cluster are located in much farther 
distance from the populations of first cluster/ 
group. Therefore we can suggest that different 
populations have gained different genetic 
variations which may be used in populations 
adaptations as also discussed in cytogenetic 
section before. 

Therefore, the present study shows 
cytogenetic and molecular diversity available 
in the C. arvense populations which is possibly 
related to their geographical distribution and 
adaptive nature of these changes, however, 
whether these genetic diversity bring about 
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morphological changes in the populations 
studied and if we can recognize any intraspecific 
taxa, need further morphological investigations 
which is under way. 
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